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ABSTRACT 



A method for predicting swell from tropical cyclones 
using a spectral wave model (TYWAVES) was tested. The model 
was applied to predicting swell propagating from three 
typhoons in the Western North Pacific through gaps in the 
Ryukyu Islands into a region of the East China Sea. The 
model involves a source region concept which considers 
only the swell emanating from regions of peak energy in 
moving typhoons. For three representative typhoons, 
predicted heights were not significantly different from 
the observed heights. The time of occurrence of the pre- 
dicted peak height agreed well with observational values 
for the swell from two typhoons, but lagged by 6-12 hours 
for the third. 

The dominant swell period and direction predicted 
by the model were not verifiable by data available for 
this study. 

Shoaling and refraction effects were considered in the 
prediction, in a sim.plified way, but attenuation was 
ignored even for the passage of energy through the Ryukyu 
Islands . 
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I . INTRODUCTION 



A. OBJECTIVE OF THE STUDY 

A serious problem for naval activities, ocean industry, 
shore protection and fisheries in the Korean south coastal 
area is the lack of adequate estimation of the waves re- 
sulting from typhoons. The objective of this study is to 
tieke a partial test of the TYWAVES for forecasting swell 
from tropical cyclones which arrive at a single observa- 
tion site near Cheju-do, Korea (33.2°N 126. 6°E). 

Forecasting based on wave fields predicted in a 
typhoon area by the TYWAVES [13] developed by NEPRF 
(Naval Environmental Prediction Research Facility) are to 
be verified against swell observations made near Cheju-do. 

B. TYPHOONS IN THE WESTERN NORTH PACIFIC 

During 1979, the Western North Pacific experienced 
28 tropical cyclones. Table I, from "1979 Annual Typhoon 
Report" [2] , shows the significant tropical cyclones for 
that year. Table II from [2] shows the monthly distribu- 
tion of tropical cyclones in 1979 and other statistics. 

Most typhoons occur in the summer season and their 
tracks can be classified as one of three typical tracks. 

The first type is that passing south of Taiwan toward the 
west, the second is that crossing over Korea through the 
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Ryukyu Islands , and the third is that passing east of 
Ryukyu Islands. 

I selected one typhoon of each type; one each in July, 
August, and September in 1979, They are Typhoons Hope, 
Irving and Owen shown in Figures 1, 2, 3 and 4 from [2]. 
tropical cyclone "best track" information is shown in 
Tables III, IV, and V for Typhoon Hope, Irving and Owen, 
respectively. 

C. TYWAVES MODEL (THE TYPHOON WAVES PROGRAM) 

The SOWM (Spectral Ocean Wave Model) run at FNOC 
(Fleet Numerical Oceanography Center) utilizes a coarse 
operating grid system (100~190NM) which does not allow 
sufficient resolution to describe adequately the resultant 
sea state in typhoon areas. Thus, TYWAVES, an improved 
model for typhoons with a locally finer grid, was developed 
by NEPRF. FNOC judged that TYWAVES is worthy of evalua- 
tion as a possible operational typhoon sea state model 
[14]. 

The TYWAVES is intended to produce fields of signifi- 
cant wave height and spectral wave properties on a mesh 
size consistent with the scales of tropical cyclones and 
is designed, primarily, for the application to the 
western North Pacific, 

The detailed outputs of TYWAVES are in the form of 
fields, at each of 12 points, of spectral energy 
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components, significant wave heights, maximum wave periods, 
and the predominant wave directions, for 00, 12, 24, 48 and 
72 hours, where 00 hours is the time of the first typhoon 
warning issued. An example of those outputs for Typhoon 
Owen is shown in Appendix A. 

Using these spectral energy components at selected 
source points, I was able to make predictions based on the 
propagation of these components as swell into the region, 
south of Cheju-do, Korea (see Map I, page 21). 

More details of TYWAVES are described in Ref. 11. 

D. FORECASTING AND VERIFICATION 

Extremely high sea states are known to be generated in 
the quadrant to the right of the direction of movement of 
typhoons. The wave generation in that quadrant of the 
typhoon derived in TYV7AVES from a spectral model utilizing 
the Pierson-Moskowitz (1964) spectrum. The model describes 
the spectroangular components of the waves present at a 
n\mber of grid points in the region of strong winds. Each 
spectral component of interest is then permitted to 
propagate at its appropriate group velocity to the fore- 
cast site. The method is applied to three western North 
Pacific typhoons in 1979 and the forecast products are 
compared with the observed swell data from the south 
coastal area of Cheju-do, Korea. 
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1979 SIGNIFICANT TROPICAL CYCLONE STATISTICS 
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FORMATION ALERTS 



23 of the 27 (85%) Formation Alert Events developed Into tropical cyclones. 
5 of the 28 (18%) tropical cyclones did not have a Formation Alert. 



WARNINGS 



Number of warning days; 149 

Ntnber of warning days with 2 cyclones; 38 

Number of warnincj days with 3 or more cyclones: S 
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Figure 3. 
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The dispersive arrival at the prediction site of the 
waves propagating from a single source region is repre- 
sented by a plot of period versus arrival time in Figures 
9, 10 and 11. The energy in each low frequency (T'lS sec) 
spectral component of the swell at Cheju-do is calculated 
by modifying the energy spectrum at the source for the 
refraction and shoaling. The lowest period (T = 4 sec) 
is assumed totally eliminated by attenuation. 

The total propagated energy at any given time is 
estimated by summing all the directional components. 

The predicted swell heights H 1/10 assumed to be related 

to the energy by H 1/10 =k/E , where K is constant (0.4) and 

6 2 

E is energy in a unit of 2 • kiloerg/cm and H 1/10 is in 
meters . 

By graphing the energy associated with the 5 mean 
period bands as a function of arrival time at the forecast 
site and summing, then a plot of the total propagated 
energy is obtained (see Figures 13, 15 and 17) . This 
plot is compared with the observed wave heights in 
Chapter IV-D. 

E. TERMINOLOGY 

To make clear the notation used in this study, the 
following terms and symbols are defined. 
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sea 



swell 



R 



o 



R 



T 

Hl/10 



H 1/3 



tar 

t 

L 

eg 

c 

D 

g 

a 

k 

h(d) 



Those waves found within a generating area. 
Those waves which have moved out of a 
generating area. 

The distance from the source to the "bottom 
line" where the corresponding periods feel 
"bottom" i.e., the swell travel distance 
in deep water. The depth at the "bottom 
line" is h = 

The distance from the "bottom line" to the 
forecast site. 

The mean period in a spectral band (T = 1/f) 
The average heights of the 1/10 highest 
waves. 

The average heights of the 1/3 highest 
waves (significant wave height) . 

The swell arrival time (GMT) . 

The swell travel time. 

Wavelength (L^ in deep water) . 

Group velocity (Cg^ in deep water) . 

Phase speed (Co in deep water) . 

The total swell travel distance (D = R^ + R) 

Acceleration of gravity 

Wave angular frequency (2)C/L) . 

Wave number (27C/L) . 

The ater depth. 
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a© “ The angle between wave crest and bottom contour 

at the depth, h = . 

2 

T^ ~ The total energy at grid point 4 in the generat- 
ing area. 

- The directional energy component at grid point 4 
in generating area. 

Wave spectrum (energy spectrxim) - The distribution of 
either wave height or energy with period. The 
potential energy of the sea surface is proportional 
to the mean of the square of the wave height. 

Dispersion - A process which leads to longitudinal 

spreading of the wave energy as, in deep water, 
energy in each spectral component is propagated 
at a characteristic group velocity Cg, the long 
waves having the larger group velocity. 

0 - The mean directional bands of 16 unit points 

rose. 

- The wave height before refraction in deep 
water. 
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II. PREDICTION OF TROPICAL STORM WAVES BY TYWAVES 



A. LOCATION OF THE SELECTED POINT-SOURCES 

The TYWAVES program produces a printout of the complete 
energy spectrum for 12 points around each typhoon. Each 
point lies at a chosen distance from the storm center in 
a fixed direction, no matter what the storm movement direc- 
tion is. These points are representative origins of the 
wave energy emerging from the typhoon. 

Only those points are selected which are in the region 
of maximum winds and where important energy components are 
directed toward the distant prediction site. 

The arrangement of points whose wave spectra are used 
in this study is shown below. 





+ 9 




+ 5 


+ 1 


+ 6 


+ 3 




+ 4 


+ 7 


+2 


+ 8 




+10 





The distances from the storm center to points 1, 2, 3 
and 4 are three grid-lengths (3 x 40 = 120 NM) , to 9, 10, 
11 and 12 are five grid lengths (5x40= 200 NM) , and to 
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5,6,7 and 8 are 3 x /2 grid lengths (3 x /2 x 40 = 170 NM) , 
respectively. 

Among the 12 points, only points 4, 6, 8 and 12 were 
selected as sources, because the spectral energy from only 
these 4 points can propagate to Cheju-do, Korea (33.2°N 
126. 6°E) . 

The locations of the sources for each time of interest 
for the three storms are shown below in Table VI, which 
was derived and calculated from the postanalysis "best 
track" data shown in Tables III, IV, and V [2]. The 
relationship between latitude and longitude is: 

Lat. = Long, x cos (Lat.) 

= Long. X cos { (storm center + forecast)/2} 

lat 

e.g. : 

When the storm center (best track) is 20.0^N 126. 7°E, 
point 4 is same lat (20.0°N) , but 120 NM (2° in lat) 
east of the center. 

Thus, 2 Lat = Long x cos ( ^ > 

= °Long X cos 10*^ 

° Long = 2° Lat/cos 10° = 2.12 8° E 

= 126.7 + 2.128 
= 128.8° E 

Therefore, point 4 ~ 20.0°N 128. 8°E. 
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Typhoon Centers and Selected Points Location 
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B. SEA-SPECTRA AT THE POINT SOURCES 



The description of the wave fields using the spectral 
component method developed by NEPRF was used in this study. 
Only those directional components able to propagate to the 
forecast site need be considered. 

A schematic example of the TYWAVES spectral prediction 
in a typhoon area is given below. The whole spectral table 
for the selected sources is shown in Appendix B. 



|t - s pectrum 

^ - mean 
i period bands 



T - 0 spectrum 
(80 components) 
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Table VII is one example of a period-direction spectrum 
at source point 4 in Typhoon Owen, at 06 on 25 September 
1979. 

The mean period is 13 sec with the maximum energy 86 
(2^- kiloeng/cm^) in the dominant period. The dominant 
vave direction is SE, but the energy from that direction 
will not arrive at the forecast site, only the SSE component 
will arrive at the site (as seen in Figure 5) . 
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C. SELECTION OF THE DIRECTIONAL ENERGY AT POINT SOURCES 
A constant energy was asstuned between parallel rays 
separated by width 80-160 NM (2-4 grid length) , which is 
40-80 NM right at the grid point depending upon the typhoon 
size. This is why the energy in each component is assumed 
constant along the deep water propagation path. As example 
of possible refraction diagrams under the above considera- 
tion and bottom contours in Map 1 is shown in Figure 5. 

Only the SSE component was considered to propagate to 
the forecast site. The SE component would refract toward 
China, west of Cheju-do, and the S component would pass to 
the east of Cheju-do as seen in Figure 5. 

All the directional components were chosen in the same 
way, and are underlined as shown in Table 7. 




Figure 5. Possible 
refraction diagram at 
point 4 at 06 GMT 
September 25, 1979. 

The curves around 
typhoon center are 
contours of wave height 
(H 1/3 in feet) . The 
detailed bottom contours 
are shown on Map 1. 
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TABLE 7 



Period-Directional Spectrum at Point 4 at 06 GMT September 
25, 1979. The energy units are 2^* kiloerg/cm2. 
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D. TOTAL AND DIRECTIONAL ENERGIES AT POINT SOURCES 

Figures 6, 7 and 8 show the energy components at each 
selected grid point in the typhoon area. Each line is 
labeled by selected grid point (source point) . T^ indicates 
total energy given by TYWAVES at point n and D^ represents 
the sum over all periods of the energies with proper 
direction to reach Cheju-do. 
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For example, v/here n = 12, 6 is SE, then 

Di2(9=SE) = E^ (T=4, 6^=SE) + E^_q (T=10, 6q=SE) 

+ ^13 (T=13, 0q=SE) + E^g (T=16, 0^=SE) 

+ E^g (T=19, e^=SE) . 

This represents the total energy at point 12 directed 
from SE toward Cheju-do. 

The energy is in practice converted into height H 1/10 
by making H 1/10 = 0.4 /e, which is described in Chapter I.D. 
Some special features of the wave fields in each source are 
now discussed. 

1 . Typhoon Hope 

The maximiam sustained wind speed was over 100 kts 
for two days from OOcMTJuly 31. The maximum wave height (H 1/10) 
recorded was 12.80 m. The directional energy components from 
grid point 12 were ignored from 00 GMT July 31 to 00 GMT 
August 01, because they were too small. 

2 . Typhoon Irving 

Only the propagation of energy components until 
00 GMT August 16 was considered, because the forecast site 
is inside of typhoon area after that time. Therefore, after 
00 GMT August 16, the propagated wave energy should be 
combined with local energy at prediction site. 
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Among the directional components, only those from 
grid points 4 and 12 were propagated to the prediction site 
during the typhoon period, because the energies from other 
sources were refracted away from Cheju-do. 

3 . Typhoon Owen 

After 06 GMT September 29, the typhoon passed east 
of southern Japan. No further propagation to Cheju-do 
occurred. 

In most cases, only one directional energy component was 
selected at each grid point. But both components (SE and 
SSE) were able to reach the forecast site on 18 GMT September 
27, from grid point 4, because of the width of their ray 
boundaries at distances near the source. 
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III. PREDICTION OF SWELL CHARACTERISTICS 



The periods of high waves due to typhoons are usually 
in the range of 7-19 sec. They travel in deep water with 
the group corresponding velocities 10-30 kts. It normally 
necessitates 1-4 days for conspicuous swells generated in 
a typhoon in the far south seas to reach Cheju-do, Korea. 

The group velocities of periods 7 and 10 sec are 
constant along their paths since they remain in "deep 
water." But for periods 13, 16 and 19 sec, the group 
velocities vary because of depth variations. Thus, I 
considered two distance components: in deep water and 

in shallow water R for periods 13, 16 and 19 sec. 

The travel distances between grid points and observa- 
tion points are calculated with the assumption of plane 
geometry. Then, group velocities, travel distances and 
travel times are derived for each period in this section. 

A. THE SETTING OF THE FORECAST SITE AND ITS RELATION TO 

THE WAVE SOURCES 

The passes of the Ryukyu Islands act as windows between 
the energy sources and the forecast site, Cheju-do, Korea. 
Since the largest land length for blockage of energy is 
only about 40 NM, this length is not sufficient to inter- 
rupt totally the energy propagation. Thus, as a further 
simplification, I neglected the effect of the Ryukyu Islands 
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against the energy propagation. The window is assumed 
sufficient for total propagation. 



B. GROUP VELOCITY 

Some simplications are used in assigning group velocity 
values along the route from typhoons to the observation 
site. 

From linear wave theory group velocity, C is given by 

U 

c - I ^ 2kh n 
^g 2 ^ sinh 2kh^ 

(all the symbols are defined below) 

and phase speed c is given by c = ^ tanh kh. In deep water 

this is well approximated by c^ = g/o = 3.03 T (ktc) for 

T given in seconds, and in shallow by c = But in the 

general case c varies with both the depth of water and the 

wavelength . 

The classifications "deep” and "shallow” are given in 
terms of "relative depth h/L^, described below. Since 
there is no "shallow water" between the wave sources and 
Cheju-do, Korea, (corresponding to the periods used by 
TYWAVES) , I will use only the general equation and the 
deep water approximation. In these equations 
Cg = group velocity 
c = phase speed (c^ in deep water) 
n = the fraction of energy propagated at phase 
speed 
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g = gravity acceleration 
h = water depth 

a = v/ave angular frequency (2T/T = 2kf, where T is 
period) 

k = local wave number (2TT/L, where L is wave length) 
f = wave frequency (1/T) 

It is important to note that c and group velocity, Cg 
must be found as functions of both h and a. A straight- 
forward method is illustrated below, where the dependence on 
kh is replaced by the deep water wavelength = 27T q/n^. 

The local wavelength is then found from L/L^, a function 
of the relative depth h/L^, 

Since C/C^ = b/L^, the phase speed c can be found from 
h and a. 

To estimate the group velocity by definition, C = 

y 

4 c [1 + — r^PVir] = nc where, n takes on the following 
2 sinh kh 

values for corresponding h/L^: 

general case (j ^ ' l>n>-^ , 

deep water case (h/L^ > ^) ; n = j . 

In order to simplify calculations, I used n = 3/4 in 
the waters between the Ryukyu Islands and Cheju-do 
where, for the longer periods, the relative depth range 

4 >h/L >-^ . With this simplication , the corresponding 
2 “ o~20 

group velocities become the following: 
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mm 







T = 13 sec; = 29.4 kts 

T = 16 sec; C = 33.6 kts 

y 

T = 19 sec; = 35.8 kts 

The influence on the calculation of C of the use of 

g 

these representative constant values is discussed in 
Chapter IV. D. 

For deep water 

= 3.03T (kts for T in sec) 

C = nc = ^ C = 1. 515T (kts for T in sec) 

g© ^ ° 

Therefore, group velicity in knots and depth at which 
Cg replaces are shown below: 

T (sec) 7 10 13 16 19 

C /C (kts) 10.6/- 15.2/- 18.7/29.6 24.3/33.6 28.8/35. C 

^o ^ 

Lq 

h(- 2- in ft) 126 256 433 656 924 

C. TRAVEL DISTANCE AND ARRIVAL TIME 

The travel distance between grid points in the storm and 
the observation point 33.2°N 126. 6°E are calculated with the 
assumptions of plane geometry. Consider this example. 

On 12 GMT July 30, the grid point 4 was at 17. 4^^133 . 

Thus, the north-south component is given by 
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y - [lat (site)° - lat (grid point) °] x 60 NM 
= [33.2 - 17.4] X 60 
= 948 NM 

and the east-west component is approximated by 

X = [long (grid point) - long (site)] cos (mean lat)] 

X 60 NM 

= [133.9 - 126.6] cos (17.4) x 60 NM 
= 418 NM 

Thus, the travel distance, D, is given by 
2 2 ^ 

D = (x + y ) ^ and for grid point 4, 

D = (418^ + 948^)^ = 1036 NM 
The travel time, t, is given by 

t = ^ + tP- , where R + R = D 

C ^ O 

R^ = the deep water distance 
R = the remaining distance 

Therefore, the arrival time, tar of each period band is 
given by 

tar (T) = t^ + t, where t^ = leaving time from 

typhoon area 
t = travel time 

The Figures 9, 10 and 11 show the swell arrival time 
for each of the three typhoons. Each curve indicates the 
arrival time of energy from a source at a specified leaving 
time labeled with various periods (from 7 to 19 sec.). 
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Appendix C shows the arrival time calculations for 



each typhoon. 
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Arrival time (GMT) 
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Arrival time (GMT) 
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Arrival time (GMT) 



IV . THE ENERGY PROPAGATION PROCEDURE 



A. PROPAGATION OF SPECTRAL ENERGY COMPONENTS 

The spectral forecast permits the forecaster to deal 
with only that range or periods which have important energy. 
Each spectral component is tracked with its respective group 
velocity , only those directions being chosen for which waves 
can reach Cheju-do. When the typhoon is moving with a 
velocity component toward Cheju-do, only spectral components 
with a group velocity greater than the movement of the 
typhoon are considered (as in Typhoon Irving) . 

Wave energy generated in a relatively small region at 
all frequencies will spread over a much larger region as it 
propagates outward from its source, and the wave character- 
istics change in such a way as to become more "swell like." 
There are essentially three processes which contribute to 
this change in wave characteristics: dispersion and angular 

spreading, which are modified by shoaling and refraction, 
and attenuation. 

In the procedure followed here, I have ignored attenua- 
tion for waves of long period. Dispersion and angular 
spreading are accounted for by simply following components 
to Cheju-do with appropriate shoaling and refraction 
factors . 
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processes 



For ease of calculation shoaling and refraction 
are simplified as described below. As seen on Map 1, a 
slightly curved contour of 100 fathoms (590 ft) connects 
the northern tip of Taiwan to the south-western tip of 
Japan. The slope from 590 ft bottom line to 413 ft (70 
fathoms) is very steep. But most bottom topography along 
the path to Cheju-do, north of the 413 ft bottom contour 
IS almost flat with the depth of 45 fathoms (266 ft) up to 
the forecast site (33.2°N 126. 6°E). 

Therefore, I have approximated to underwater topography 
by assuming only two water depths, a deep water and a 
shallow water (intermediate water depth) region with an 
abrupt jump between them. 



To compute the wave characteristics at a shoal water 
site, the shoaling and refraction are considered using 
the values of and n from Chapter III.B. 

The energy of component (T,9) in the typhoon area 
(deep water) is transformed after shoaling and refraction 
to its value at the forecast site according to 
E(i, j) = E^ (i, j)-kj (i, j) kj (i, j) 
where E (i, j) is the energy of the component of 
period T^ and which had the direction 9^ in the 
generating area, 

and kg (i, j) , k^ (i, j) are the respective 
shoaling and refraction factors of those components. 
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Details of calculation of k and k are given in the 

s r 

following sections. The refraction and shoaling calcula- 
tions for three typhoons are shown in Appendix D. The 



following is 


a sample 


calculation 


on 12 GMT 


July 30 


for 


Typhoon Hope 


at source 


point 4. 








T 


7 


10 


13 


16 


19 


k 

s 


1 


1 


0.92 


0. 89 


0.90 


k 

r 


1 


1 


1 


0.989 


0.972 


k 

s r 


1 


1 


0.8464 


0.7748 


0.7653 


E 

o 


7 


19 


11 


9 


14 


E 


7 


19 


9.3 


7.0 


10.7 



where k is derived from plate C-1 [6] 
with h = 40 fathoms = 236 ft 

is derived from Figure 2-19 [5] 
with h = 413 ft and = 40°. 

B. SHOALING AND REFRACTION FACTOR 
1 . Shoaling Factor 

The shoaling factor, = H/H^ was derived from the 
plate C-1 [6] with depth and period of the appropriate 
spectral component. The shoal site depth near Cheju-do 
(33.2°N 126. 6°E) is 40 fathoms, and the relative depth h/L^ 
and shoaling factor of each period are shown below. 
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Figure 2-19 Changes m Wove Direcfion and Height Due to Refraction on Slopes with 
Straight , Parallel Depth Contours 



T 


17 


10 


13 


16 


19 


h/L 

o 


0.94 


0.46 


0.24 


0.18 


0.13 


k 

s 


1 


1 


0.92 


0.89 


0.90 



2 . Refraction Factor 

Generally, two basic techniques of refraction 

analysis are available — graphical and numerical. Several 

graphical procedures are available, but all methods of 

refraction analyses are based upon Snell's law. Refraction 

may be treated analytically in any region with straight 

and parallel contours, by using Snell's law directly: 

sin a = (— — ) sin a , where a is the angle between the wave 
^o ° 

crest and the bottom contour, and is the angle between 

the deep water wave crest and the bottom contour. 

Figure 2-19 [5] shows the relationships among 
a, a^, period, depth and refraction factor in graphical 
form. I derived the refraction factor from using this 
graph. Figure 2-19, the bottom contour and period. Also, 

I assumed that the refraction factor is 1 if the angle 
(between crest and bottom contour) is less than 10*^ and 
that the refraction occurs only one time at the depth of 
maximum bottom gradient where waves are refracted because 
of the wide flat-bottomed portions over most of the inter- 
mediate water propagation path as seen on Map 1. 
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with these considerations, I derived the refraction 
factor and predicted the shoaling energy. Those procedures 
are shown in Appendix D. 

C. SHOALING AND REFRACTION OF THE SPECTRAL COMPONENTS 

Each deep water wave spectral component derived from 
the TYWAVES model moves with its respective group velocity. 
I considered the 80 energy components to behave as mono- 
chromatic component waves. To assess the shoaling and 
refraction of each to the shoal-water site I assumed that 
the wave power transmitted between a given pair of 
orthogonals remains constant at all depths (this means no 
frictional losses, diffraction or scattering, and also 
implies that a steady state exists) . 

With these assumptions the wave power P is given by 
P = EC b=E^C b^ 

g o o 

C b^ 

Thus, E = == ■'s = I 

g 

where E = the average wave energy per unit area 
of sea surface for waves transformed by shoaling and 
refraction. 

P = water density 
g = acceleration of gravity 
H = wave height of transformed waves 
Cg ,^g ~ group velocity 
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b^,b = orthogonal separations 
= shoaling factor 
= refraction factor 

The energy in each spectral component of the swell at 
the observation site was calculated by modifying the energy 
spectrum at the selected sources for the effects of shoaling 
and refraction according to equation E = E^*k^*k^. 

The energy associated with the various components as a 
function of time of arrival at observation site, as seen in 
Figures 9, 10 and 11. 

The total energy in the swell at any given arrival time 
is estimated by summing all the shoaling components at that 
time . 

In summary, each component in typhoon area's T - 6 
spectrum is shoaled and refracted using the k^ and k^ 
values appropriate to it to find the energy at the observa- 
tion site. 

1. Shoaling Process (Computation) 

See Appendix D. 

2 . Shoaling Energy Components From Each Source Versus 
Arrival Time 

As seen in Figures 12, 14 and 16, the energy of 
the components from all sources is shown as a function of 
its arrival time at Cheju. See Appendix D. 
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3. The Predicted Swell Waves at Prediction Site 



As shown in Figures 13, 15 and 17, the predicted 
swell waves are the sum of all transformed components at 
given time at the forecast site. 

For Typhoon Irving (see Figure 15) beginning 
12 GMT August 16, the forecast site is already inside 
typhoon area. Therefore, there is no prediction done 
after that time. These predictions are discussed in 
Chapter IV. D. 

D. THE OBSERVED DATA 

Observations of wave conditions for verification of 
the swell forecasts were obtained from the sources, the 
National Oceanographic Data Center (NODC) and the Republic 
of Korea (ROK) Navy. All listed heights are based on 
visual estimates. 

Following Table VIII shows the visual observation of 
swell heights (H 1/10) for each typhoon. The data from 
NODC are sparse and often far from the forecast site, 
but they appear to be samples from the same set as those 
of the ROK Navy. The data from ROK Navy visual observa- 
tions were made at 33.2°N, 126. 5°E close to my point of 
interest (33.2°N, 126. 6°E). So I considered this point 
exactly the same as my forecast site. 
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Arrival time (GMT) 



TABLE 8 



The Visual Observation Data by NODC and ROK Navy 
1. For Typhoon Hope 



DTG (GMT) 


Location (°N-°E) 


H 1/10 (meters) 


Source 


79080100 


33.6-125.2 


2.5 


NODC 


0104 


32.3-126.0 


2.5 


NODC 


0112 


31.2-126.7 


4.0 


NODC 


0118 


27.7-124.6 


6.5 


NODC 


0121 


27.8-125.4 


6.5 


NODC 


0200 


29.1-128.0 


6.5 


NODC 


0206 


28.0-129.0 


5.5 


NODC 


0218 


26.9-125.8 


5.0 


NODC 


0300 


28.7-128.6 


3.0 


NODC 


0306 


31.4-126.2 


3.0 


NODC 


0312 


32.3-125.9 


3.5 


NODC 


080100 


33.2-126.5 


3.0 


ROK Navy 


0104 


33.2-126.5 


3.5 


ROK Navy 


0108 


33.2-126.5 


3.7 


ROK Navy 


0112 


33.2-126.5 


4.0 


ROK Navy 


0116 


33.2-126.5 


4.5 


ROK Navy 


0120 


33.2-126.5 


4.5 


ROK Navy 


0208 


33.2-126.5 


5.0 


ROK Navy 


0212 


33.2-126.5 


4.5 


ROK Navy 


0216 


33.2-126.5 


5.5 


ROK Navy 


0220 


33.2-126.5 


4.5 


ROK Navy 


0308 


33.2-126.5 


3.5 


ROK Navy 


0312 


33.2-126.5 


3.5 


ROK Navy 


0316 


33.2-126.5 


3.0 


ROK Navy 


0320 


33.2-126.5 


5.0 


ROK Navy 


0408 


33.2-126.5 


4.0 


ROK Navy 


0412 


33.2-126.5 


3.0 


ROK Navy 


0416 


33.2-126.5 


3.0 


ROK Navy 


0420 


33.2-126.5 


2.5 


ROK Navy 
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2 . For Typhoon Irving 



DTG (GMT) 


Location (°N-°E) 


H 1/10 (meters) 


Source 


79081400 


29.10127.6 


2.5 


NODC 


1403 


28.5-127.6 


4.0 


NODC 


1406 


28.7-129.3 


3.5 


NODC 


1409 


27.7-129.7 


4.0 


NODC 


1412 


29.2-131.1 


4.5 


NODC 


1500 


30.0-125.0 


6.0 


NODC 


1503 


30.7-127.3 


4.0 


NODC 


1515 


28.3-130.3 


8.0 


NODC 


1518 


27.2-130.0 


9.0 


NODC 


1521 


26.5-128.8 


8.0 


NODC 


1600 


27.1-130.0 


6.5 


NODC 


1603 


27.3-128.5 


8.0 


NODC 


1606 


29. 6-129.2 


5.5 


NODC 


1700 


27.4-125.5 


6.0 


NODC 


79081400 


33.2-126.5 


2.5 


ROK Navy 


1408 


33.2-126.5 


2.5 


ROK Navy 


1412 


33.2-126.5 


3.0 


ROK Navy 


1416 


33.2-126.5 


3.5 


ROK Navy 


1420 


33.2-126.5 


3.5 


ROK Navy 


1508 


33.2-126.5 


4.0 


ROK Navy 


1512 


33.2-126.5 


4.5 


ROK Navy 


1516 


33.2-126.5 


5.0 


ROK Navy 


1520 


33.2-126.5 


6.0 


ROK Navy 


1608 


33.2-126.5 


8.0 


ROK Navy 


1612 


33.2-126.5 


9.0 


ROK Navy 


1616 


33.2-126.5 


8.0 


ROK Navy 


1620 


33.2-126.5 


8.0 


ROK Navy 


1708 


33.2-126.5 


8.0 


ROK Navy 


1712 


33.2-126.5 


9.0 


ROK Navy 


1716 


33.2-126.5 


8.0 


ROK Navy 


1720 


33.2-126.5 


7.0 


ROK Navy 
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3 . For Typhoon Owen 



DTG (GMT) 


Location (°N-°E) 


H 1/10 (meters) 


Source 


79092600 


31.4-128.0 


4.0 


NO DC 


2606 


33.6-129.1 


3.5 


NODC 


2615 


31.8-129.5 


3.0 


NODC 


2700 


28.1-126.9 


6.0 


NODC 


2703 


27.1-126.6 


6.5 


NODC 


2706 


32.6-125.7 


4.0 


NODC 


2709 


31.6-125.7 


5.5 


NODC 


2712 


30.7-129.8 


6.5 


NODC 


2715 


28.7-125.3 


6.0 


NODC 


2800 


32.8-128. 3 


5.5 


NODC 


2803 


30.6-125.7 


6.5 


NODC 


2806 


31.5-126.8 


6.5 


NODC 


2809 


31.6-126.6 


6.0 


NODC 


2812 


31.7-126.6 


6.5 


NODC 


2818 


34.0-128.3 


6.0 


NODC 


2900 


32.7-127.7 


6.5 


NODC 


2606 


32.1-127.3 


5.5 


NODC 


2912 


31.0-126.4 


5.5 


NODC 


2918 


34.0-129.7 


5.0 


NODC 


3000 


29.0-124.7 


4.0 


NODC 


3006 


30.8-127.7 


3.0 


NODC 


3018 


30.4-131.6 


4.0 


NODC 


79092600 


33.2-126.5 


2.5 


ROK Navy 


2608 


33.2-126.5 


4.0 


ROK Navy 


2612 


33.2-126.5 


5.0 


ROK Navy 


2616 


33.2-126.5 


5.0 


ROK Navy 


2620 


33.2-126.5 


4.5 


ROK Navy 


2708 


33.2-126.5 


5.5 


ROK Navy 


2712 


33.2-126.5 


5.0 


ROK Navy 


2716 


33.2-126.5 


5.0 


ROK Navy 


2720 


33.2-126.5 


5.5 


ROK Navy 


2808 


33.2-126.5 


4.5 


ROK Navy 


2812 


33.2-126.5 


5.0 


ROK Navy 


2816 


33.2-126.5 


5.5 


ROK Navy 


2820 


33.2-126.5 


5.5 


ROK Navy 


2908 


33.2-126.5 


4.5 


ROK Navy 


2912 


33.2-126.5 


4.0 


ROK Navy 


2916 


33.2-126.5 


4.5 


ROK Navy 


2920 


33.2-126.5 


4.5 


ROK Navy 


3008 


33.2-126.5 


5.0 


ROK Navy 


3012 


33.2-126.5 


4.0 


ROK Navy 


3016 


33.2-126.5 


2.5 


ROK Navy 


3020 


33.2-126.5 


3.0 


ROK Navy 
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E. COMPARISON OF SWELL PREDICTIONS AND OBSERVATIONS 

The predicted and observed heights are plotted in 
Figures 13, 15 and 17, and are good agreement. Detailed 
comparisons for each typhoon follow: 

1 . Typhoon Hope 

The first predicted swell arrival times from each 
of three sources, and the NODC and ROK Navy observations 
are almost the same with similar energies (see Figure 13) . 

The peak energy arrival times are also nearly the 
same, but the peak predicted energy is slightly higher than 
the observations. This is reasonable agreement since it is 
estimated that the visual observations can be error. Verploegh 
1961 [18] estimated the average observational error for a 
visual observation of wave height varies from 1 ft at 5 ft 
wave heights to 3 ft at 18 ft wave heights. 

2 . Typhoon Irving 

As seen in Figure 15, the predicted swell arrival time 
nearly matches the time shown by both sets of observations. 
However, the predicted energy of the rise is slightly higher 
than the observational values. The peak energy (around the 
time of 20 GMT August 15) is also the same as the observa- 
tional peak. 

After 00 GMT August 16, the prediction ' site is within 
the typhoon's wind circulation area. Most wave heights (H 1/10) 
are 8m and the highest observation is 10m. These waves 
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dominate the entire wave field and no swell can be dis- 
tinguished later than 00 GMT August 16. 

3 . Typhoon Owen 

In Typhoon Owen (see Figure 17) , the swell arrival 
time indicated by observations is earlier than the predicted 
time by about 6-12 hours. The predicted peak energy lies 
close to the observational values. 

In fact, the NODC values are higher than those of 
the prediction, and ROK Navy observations are lower. 

Except near the peak both observation sets show similar 
time variations. Therefore, the observed values are con- 
sidered consistent and probably accurate. 

F. ERROR SOURCES 

I have made several assumptions in this study in order 
to simplify calculations. The most serious error sources 
involve assumptions about the windows in the Ryukyu Islands, 
group velocity in shallow water, and simple bottom contours 
in shallow water. 

There are also differences between the predictions (which 
include waves dependent on the local weather conditions) that 
make it difficult to evaluate the prediction. Regarding the 
assumptions about group velocity in Chapter III.B, I used 
n = 3/4 in shallow water for calculation of group velocity 
Cg = nc. In Typhoon Hope, which has the longest shallow 
water travel distance, miscalculation of C would have its 
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greatest effect. Yet the observed and predicted energy 
peaks are not greatly separated. This may mean that the 
approximations are realistic. 

The bottom contour assumptions seem to be reasonable 
through three tests. The local weather condition is the 
most serious factor contributing to differences between 



the observations and the predictions. The local wind 
pattern was the following at Cheju-do, according to the 
ROK Navy data set. 



1979 072900-073100: SW 5-8 kts (preceding Hope 

swell) 

080100-080300: SW 6-10 kts (during Hope 

swell) 

1979 081200-081300: S 6-8 kts (preceding Irving 

swell) 

081400-081700: S 30-60 kts (during Irving 

swell) 

092400-092600: E 5-8 kts (preceding Owen 

swell) 

092600-093000: E 20-30 kts (during Owen 

swell) 

Therefore, before the swell arrivals, the local wave heights 
were considered less than 1~1.5 m. During Irving and Owen 
there were important local sea contributions to observed 
height. 

Lastly, the original error sources, i.e., those involving 
input parameters for the TYWAVES model, the location of the 
typhoon center, the winds in the source region, and typhoon 
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Those errors are 



size, etc., are ignored in this study, 
discussed in Refs. [2] and [17] . 
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V. CONCLUSIONS 



From Figures 13, 15 and 17, I can conclude that most 
predicted swell heights are lower than those of the obser- 
vations (combined sea and swell) . 

The times of occurrence of the predicted peak heights 
agreed reasonably well with those of observations for the 
swell from each typhoon. 

These results on the basis of this limited test suggest 
that TYWAVES predict satisfactorily those swells in the East 
China Sea which originate in tropical cyclones in the 
western North Pacific. 

Computer aided predictions may improve the quality 
of the forecasts by reducing the need for simplifying 
approximations, as in the case of the treatment in this 
thesis of the shoaling and refraction processes. Such 
predictions would also provide a larger base for assessing 
the accuracy of the method. Additional verification is 
required to draw more specific conclusions. 
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APPENDIX A 



Outputs of TYWAVES For Typhoon Owen 
at 12 GMT September 26 

The following tables and figures show the computer anal- 
ysis for Typhoon Owen by TYWAVES. 

a. Period-directional spectrum at each of 12 sources. 
b-1. The significant wave height (H 1/3 in feet) distri- 
bution in typhoon area. "0" indicates the land 
area and the distance between the grid points is 
4 0 NM. 

b-2. The maximum wave period distribution in typhoon 
area. "-i" indicates the land area. 
b-3. The maximum wind wave directions distribution. 

The direction indicated with 16 unit point rose, from 
North (1) to NNE (16) with CCW direction, 
c. NEPRF Typhoon Wave Program Analysis for Typhoon Owen 
during the typhoon period from the time of the first 
typhoon warning to 72 hours later, based on post- 
analysis data. 
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TOTAL 272 J 9 6 1 52 51 TOTAL 172 0 13 60 58 22 17 TOTAL 3.«0 

periods 4 7 Ij 13 16 19 PERIODS 4 7 10 13 16 19 PERIODS 
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The significant wave height (H 1/3 in feet) distribution in typhoon area. "O" indicates 
the land area and the distance between the grid points is 40 NM. 
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The maximum wave period distribution in typhoon area. "-1" indicates the land area. 
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The maximum wind wave directions distribution. The direction indicates with 16 unit 
point rose, from North (1) to NNE (16) with CCW direction. 



NEPRF Typhoon Wave Program Analysis for Typhoon Owen during the 
typhoon period from the time of the first typhoon warning to 
72 hours later, based on post-analysis data. 
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APPENDIX B 



Sea Spectra at Source Region for Three Typhoons and Significant 

Wave Distribution 

The following tables and figures show the period-direction 
spectrum at selected source grid point and NEPRF typhoon wave 
program analysis for wave height (H 1/3 in feet) distribution 
around the typhoon center, respectively. 

a-1. NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected grid point) for Typhoon Hope. 

a-2. NEPRF Typhoon Wave Program Analysis for Typhoon Hope. 

b-1. NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected source grid point) for Typhoon 
Irving. 

b-2. NEPRF Typhoon Wave Program Analysis for Typhoon 
Irving. 

c-1. NEPRF Typhoon Wave Model (period-direction spectrum 
at each selected source grid point) for Typhoon 
Owen. 

c-2. NEPRF Typhoon Wave Program Analysis for Typhoon 
Owen. 
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NEPRF Typhoon Wave Model 
(Hope ) 
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NEPRF Typhoon Wave Model 
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APPENDIX C 



Propagation Energy Arrival Time 

The following tables show the components from each 
source grid point, travel time and arrival time (GMT) . 

and R are determined from Map 1 (U.S. Navigation Chart 
No. 94027, Scale 1: 927.700 at lat 32^15') from the critical 
depth of each period. 



89 



15.2 









Cgo/Cgs 


10.6 








T 


7 


12 Jul 


30 


PT4 


Ro/R 


1036 


(GMT) 




17.4-133.9 


t 


97.74 








tar 


031344 






PT12 


PD/R 


1071 






17.4-135.3 


t 


101.03 








tar 


031702 


00 Jul 


31 


PT4 


Ro/R 


919 






18.6-131.5 


t 


86.70 








tar 


031443 






PT12 


Ro/R 


946 






18.6-132.9 


t 


89.29 








tar 


031717 


12 Jul 


31 


PT4 


Ro/R 


822 






19.6-128.3 


t 


77.55 








tar 


031731 


00 Aug 


01 


PT4 


Ro/R 


780 






20.6-125.3 


t 


73.56 








tar 


040134 






PT12 


Rd/R 


756 






20.6-126.8 


t 


71.32 








tar 


032319 


12 Aug 


01 


PT4 


Ro/R 


744 






21.5-122.2 


t 


70.19 








tar 


041011 






PT8 


Ro/R 


858 






19.5-122.2 


t 


80.94 








tar 


042059 


12 Aug 


01 


PT12 


Rd/R 


720 






21.5-123.7 


t 


67.92 








tar 


040755 


00 Aug 


13 


PT4 


Cgo/Cgs 


10.6/- 






20.0-128.8 


Ro/R 


802/- 








t 


75.66 








tar 


160340 






PT12 


Rd/R 


817/- 






20.0-130.2 


t 


77.08 








tar 


160505 



18.7/29.4 24.3/33.6 28.8/35.8 



10 

1036 

68.16 

0202809 


13 

896/140 

52.67 

011640 


16 

873/163 

40.78 

010447 


19 

866/170 

34.82 

312249 


1071 

70.46 

021028 


941/130 

54.74 

011844 


876/175 

42.08 

010605 


884/187 

35.91 

312355 


919 

60.46 

021228 


749/170 

45.83 

012150 


719/200 

34.94 

011056 


709/210 

30.49 

010629 


946 

62.24 

021414 


816/130 

48.06 

020004 


781/165 

37.05 

011303 


770/176 

31.66 

010740 


822 

54.08 

021847 


610/212 

39.83 

020350 


555/267 

30.79 

011847 


550/272 

26.70 

011442 


780 

51.32 

030319 


410/370 

34.52 

021031 


370/440 

27.43 

020326 


365/415 

24.26 

020014 


756 

49.74 

030144 


501/255 

35.46 

021128 


476/280 

27.92 

020355 


471/285 

24.31 

020019 


744 

48.95 

031257 


324/420 

31.62 

021937 


284/460 

25.38 

021323 


279/465 

22.68 

021041 


858 

56.45 

032027 


446/412 

37.86 

030152 


408/450 

30.18 

021811 


406/452 

26.73 

021444 


720 

47.37 

031122 


350/370 

31.30 

021918 


305/415 

24.90 

021254 


300/420 

22.15 

021009 


15.2/- 

802/- 

6.76 

150446 


13.7/29.4 

562/240 

38.22 

141413 


24.3/33.6 

550/252 

30.13 

140608 


23.2/35.8 

546/256 

26.11 

140207 


817/- 

53.75 

150545 


617/200 

39.80 

141548 


607/210 

31.23 

140714 


600/217 

26.89 

140253 
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12 Aug 13 


PT4 

22.0-138.2 


Rd/R 

t 

tar 


678/- 

63.96 

160358 


678/- 

44.61 

150837 


433/245 

31.49 

141929 


423/255 

25.00 

141300 


413/265 

21.74 

140945 




PT12 

22.0-129.6 


Ro/R 

t 

tar 


692/- 

65.28 

160519 


692/- 

45.53 

150932 


467/225 

32.63 

142038 


452/240 

25.74 

141345 


442/250 

22.33 

141020 


00 Ai:^ 14 


PT4 

23.5-127.2 


Ro/R 

t 

tar 


583/- 

55.00 

160700 


583/- 

38.36 

151422 


323/260 

26.12 

150207 


303/280 

20.80 

142048 


300/283 

18.32 

141819 




PT12 

23.5-128.6 


Ro/R 

t 

tar 


592/- 

55.85 

160751 


592/- 

38.95 

151457 


362/230 

27.18 

150311 


352/240 

21.63 

142138 


350/242 

18.91 

141855 


12 Aug 14 


PT4 

24.6-126.7 


Ro/R 

t 

tar 


516/- 

48.68 

161241 


516/- 

33.95 

152157 


239/277 

22.20 

151012 


196/320 

17.59 

150535 


191/325 

15.71 

150343 


12 Aug 14 


PT12 

24.6-128.2 


Ro/R 

t 

tar 


523/- 

49.34 

161320 


523/- 

34.41 

152224 


288/235 

23.39 

151123 


278/245 

18.73 

150614 


263/260 

16.39 

150424 


00 Aug 15 


PT4 

25.9-126.5 


Ro/R 

t 

tar 


438/- 

41.32 

161719 


438/- 

28.82 

160449 


155/293 

18.25 

151815 


108/330 

14.27 

151416 


103/335 

12.93 

151256 


00 Aug 15 


PT12 

25.9-128.0 


Ro/R 

t 

tar 


444/- 

41.89 

161753 


444/- 

29.10 

160513 


214/230 

19.27 

151916 


204/240 

15.53 

151532 


199/245 

13.75 

151345 


12 Aug 15 


PT4 

27.5-126.0 


Ro/R 

t 

tar 


343/- 

32.36 

162022 


343/- 

22.57 

161034 


13/330 

11.92 

152355 


-/343 

10.21 

152213 


-/343 

9.58 

152135 




PT8 

25.5-125.9 


Ro/R 

t 

tar 


464/- 

43.73 

170744 


464/- 

30.53 

161832 


144/320 

18.58 

160635 


119/345 

15.16 

160310 


109/355 

13.70 

160142 




PT12 

27.5-127.5 


RoA 

t 

tar 


345/- 

32.54 

162033 


345/- 

22.70 

161042 


105/240 

13.78 

160147 


95/250 

11.35 

152321 


85/260 

10.21 

152213 


00 Aug 16 


PT4 

29.6-126.0 


Ro/R 

t 

tar 


218/- 

20.57 

162034 


218/- 

14.34 

161421 


-/218 

7.41 

160725 


-/218 

6.49 

160629 


-/218 

6.09 

160605 




PT8 

27.6-126. 


Ro/R 

t 

tar 


337/- 

31.79 

170747 


337/- 

22.17 

162210 


20/317 

11.85 

161151 


-/317 

10.03 

161002 


-/317 

9.41 

160925 
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PT12 

29.6-127.5 


Po/R 

t 

tar 


221/- 

20.85 

162051 


221/- 

14.54 

161432 


-/221 

7.52 

160731 


-/221 

6.58 

160635 


-/221 

6.17 

160610 


06 Sep 25 


PT4 

21.3-132.4 


Cgo/Cgs 10.6 
Rd/R 784 

t 73.96 

tar 280758 


15.2 

784 

51.58 

270935 


18.0/29.4 

654/730 

39.39 

262123 


24.3/33.6 

621/163 

30.44 

261225 


28.8/35.8 

611/173 

26.05 

260803 




PT12 

21.3-133.9 


Rd/R 

t 

tar 


822 

77.55 

281133 


822 

54.08 

271205 


702/120 

41.62 

262337 


662/160 

32.00 

261400 


658/164 

27.43 

260926 


18 Sep 25 


PT4 

22.6-131.7 


Ro/R 

t 

tar 


696 

65.66 

281140 


696 

45.73 

271547 


563/133 

74.63 

270438 


513/183 

26.56 

262034 


508/188 

22.89 

261653 




PT12 

22.6-133.1 


Rd/R 

t 

tar 


731 

68.96 

281458 


731 

48.09 

271805 


610/121 

36.74 

270644 


577/154 

28.32 

262219 


574/157 

24.32 

261819 


06 Sep 26 


PT4 

23.5-131.4 


Rd/R 

t 

tar 


639 

60.28 

281817 


639 

42.04 

280002 


509/130 

31.64 

271338 


469/170 

24.36 

270622 


419/220 

20.07 

270204 




PT12 

23.5-132.8 


Rd/R 

t 

tar 


675 

63.68 

282141 


675 

44.41 

280225 


555/120 

33.76 

271546 


522/153 

26.03 

270802 


515/160 

22.35 

270421 


18 Sep 26 


PT4 

24.4-131.6 


Rd/R 

t 

tar 


542 

51.13 

282108 


542 

35.66 

280539 


417/125 

26.55 

272034 


382/160 

20.48 

271429 


371/171 

17.66 

271140 




PT12 

24.4-135.3 


Rd/R 

t 

tar 


667 

62.92 

290855 


667 

43.88 

281353 


607/60 

34.5 

280430 


562/. 05 
26.26 
272016 


562/105 

22.44 

271626 


06 Sep 27 


PT4 

25.5-131.9 


Rd/R 

t 

tar 


544 

51.32 

290919 


544 

35.79 

281747 


484/60 

27.92 

280955 


410/134 

20.86 

280252 


407/137 

17.96 

272358 




PT12 

25.5-133.4 


Rd/R 

t 

tar 


591 

55.75 

291345 


591 

38.88 

282053 


526/65 

30.34 

281220 


446/145 

22.67 

280440 


443/148 

19.52 

280131 


18 Sep 27 


PT4 

26.5-132.0 


Rd/R 

t 

tar 


496 

46.79 

291648 


496 

32.63 

290238 


430/66 

25.24 

281914 


348/148 

18.73 

281244 


345/151 

16.20 

281012 
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PT6 Rd/R 

28.5- 132.1 t 

tar 

PT12 Ro/R 

26.5- 133.5 t 

tar 

06 Sep 28 PT4 Ro/R 

27.3- 132.1 t 

tar 

PT12 Rd/R 

27.3- 133.5 t 

tar 

18 Sep 28 PT4 Ro/R 

27.8-132.1 t 

tar 

06 Sep 29 PT4 Ro/R 

28.5- 132.4 t 

tar 



406 

38.30 

290818 


406 

26.71 

282043 


306/100 

19.76 

281346 


251/155 

14.94 

280857 


547 

51.60 

292136 


547 

35.99 

290559 


487/60 

28.08 

282205 


417/130 

21.03 

281502 


460 

43.40 

300124 


460 

30.26 

291216 


400.60 

23.43 

290526 


318/142 

17.31 

282319 


511 

48.21 

300612 


511 

33.62 

291537 


414/97 

25.44 

290726 


386/125 

19.61 

290136 


436 

41.13 

301108 


436 

28.68 

292241 


371/65 

22.05 

291603 


291/145 

16.29 

291017 


416 

39.25 

302115 


416 

27.37 

300922 


361/55 

21.18 

300311 


291/125 

15.70 

292142 



245/160 

13.01 

280706 

407/140 

18.04 

281202 

315/145 

14.99 

282059 

381/130 
16. 86 
282252 

289/147 

14.14 

290608 

288/128 

13.58 

291935 
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APPENDIX D 



Shoaling of the Spectral Energy Components 

For refraction factor computation at 00 GMT August 13 

from Point 4 of Typhoon Irving, the critical water depth is 

325 feet and the angle between the swell crest and the 

bottom contour of 325 ft depth is 40 . Thus, with h/gT^ 

and the Figure 2-19 [7], the refraction factors are derived 

like the follov;ing: 

T (sec) =7 10 13 16 19 

K ) = 1 1 1 0.98 0.963 

r H 

o 

In the same way, the other refraction factors are shown in 
shoaling computation, in the following tables; 

1 . Shoaling Computation 

a. Typhoon Hope 

b. Typhoon Irving 

c. Typhoon Owen 

2 . Shoaling Energy and Arrival Times 



a. 


Typhoon 


Hope 


b. 


Typhoon 


Irving 


c. 


Typhoon 


Owen 
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1. Shoaling Canputation 







T 


7 






k 

s 


1 


Typhoon Hope 






12 Jill 30 


PT4 


kr ry 


1 


((3ir) 


17.4-133.9 


(k k 
s r 


1 




h = 413 


SSE 

o 


7 




11 

O 

0 


SSE 

S 


7 




PT12 


k 


1 




17.4-135.3 


(k k 
s r 


1 




h = 413 


SSE 

o 


8 




0 

LO 

ro 

II 


SSE 

s 


8 


00 Jul 31 


PT4 


k 

r 


1 




18.6-131.5 


(k k )2 
s r 


1 




h = 413 


SSE 

o 


4 




a = 45° 
o 


SSE 

s 


4 




PT12 


k 

r 


1 




18.6-132.9 


(k k 
s r 


1 




h = 431 


SSE 

o 


3 




a = 35 
o 


SSE 


3 


12 Jul 31 


PT4 




1 




19.6-128.3 




1 




n = 413 


s 

o 


0 




-tpO 

a = 45 
o 


s 

s 


0 


00 Aug 01 


PT4 


k 

r 


1 




20.6-123.2 


(k k 
s r 


1 




h = 236 


SSW 

o 


5 




a = 60° 
o 


ssw 

s 


5 



13 

0.92 


16 

0.89 


19 

0.90 


E (2^- kilo- 
H 1/10 (in 
erg/on^) 
meters) 


1 

0.8464 


0.989 

0.7748 


0.972 

0.7653 


E 


H 1/10 


11 


9 


14 


63 


3.17 


9.3 


7.0 


10.7 


53 


2.91 


1 


1 


0.98 






0.8464 


0.7921 


0.7779 






31 


36 


42 


147 


4.85 


26.2 


28.5 


32.7 


123.4 


4.44 


1 


0.984 


0.962 






0.8464 


0.7670 


0.7496 






20 


13 


15 


69 


3.32 


16.9 


10.0 


11.2 


57.1 


3.02 


1 


1 


0.982 






0.8467 


0.7921 


0.7811 






43 


19 


19 


123 


4.44 


36.4 


15 


14.8 


106.2 


4.12 


1 


0.984 


0.962 






0.8464 


0.767 


0.75 






30 


46 


20 


121 


4.40 


25.4 


35.3 


15 


98.7 


3.97 


0.935 


0.875 


0.825 






0.7399 


0.6065 


0.5513 






30 


19 


40 


118 


4.35 


22.2 


11.5 


22.1 


83.8 


3. 66 



10 

1 

1 

1 

19 

19 

1 

1 

28 

28 

1 

1 

15 

15 

1 

1 

37 

37 

1 

1 

23 

23 

1 

1 

23 

23 
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PT12 k 1 

r 

20.6-126.8 (kk)^ 1 

s r 

h = 413 S 0 

o 

a = 35° S 0 

o s 

12 Aug 01 PT4 k^ 1 

21.5- 122.2 (kk)^ 1 

s r 

h = 236 SSW 0 
o 

a = 30° SSW 0 

o s 

PT8 k 1 

r 

19.5- 122.2 (kk)^ 1 

s r 

h = 295 SSW 1 
o 

a = 55° SSW 1 

o s 

PT12 k^ 1 

21.5- 123.7 (kk)^ 1 

s r 

h = 295 SSW 0 
o 

a = 50° SSW 0 
o s 

b. Typhoon Irving 

00 Aug 13 PT4 

(Oyn?) 20.0-128.8 

deep SSE 3 

shoal k^ 1 

h = 325 ft (kgk^) ^ 1 

a = 40° SSE 3 

o 



1 


1 


0.98 






0.8464 


0.7921 


0.7779 






71 


98 


52 


243 


6.24 


60.1 


77.6 


40.5 


198.2 


5.63 


0.99 


0.98 


0.97 






0.8296 


0.7607 


0.7621 






33 


3 


10 


58 


3.05 


27.4 


2.3 


7.6 


28.3 


2.78 


0.97 


0.93 


0.88 






0.7963 


0.6851 


0.6272 






51 


8 


13 


98 


3.96 


40.6 


5.5 


8.2 


77.3 


3.52 


0.977 


0.95 


0.915 






0.8079 


0.7149 


0.6781 






17 


49 


31 


104 


4.08 


13.7 


35.0 


21 


75.7 


3.48 



46 


21 


27 


137 


4.68m 


1 


0.98 


0.963 






0.8464 


0.7607 


0.7512 






38.9 


16.0 


20.3 


117.2 


4.33 



1 

1 

20 

20 

1 

1 

11 

11 

1 

1 

22 

22 

1 

1 

6 

6 

39 

1 

1 

39 
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IT- 







4 









PT12 



20.0-130.2 

deep SSE 5 

shoal 1 

h = 325 (k k ) ^ 1 

s r 

a = 10 SSE 5 

o 

12 Aug 13 PT4 

22.0-128.2 

deep SSE 4 

shoal k 1 

r 

h = 325 (kgk^) ^ 1 

a = 10 SSE 4 

o 

PT12 

22.0-129.6 

deep SSE 3 

shoal k 1 

r 

h = 325 (k k ) ^ 1 

s r 

a = 10 SSE 3 

o 

00 Aug 14 PT4 

23.5-125.0 

deep S 4 

shoal k 1 

r 

h = 325 (k k ) ^ 1 

s r 

a = 4.5 S 4 

o 



38 


20 


11 


100 4.00 


1 


1 


1 




0.8464 


0.7921 


0.81 




35.0 


15.8 


8.9 


89.7 3.79 



31 


25 


13 


105 4.10 


1 


1 


1 




0.8464 


0.7921 


0.81 




26.2 


19.8 


10.5 


90.5 3.81 



20 


15 


4 


67 


3.27 


1 


1 


1 






1 


0.8464 


0.7921 


0.81 




16.9 


11.9 


3.2 


59 


3.07 



22 


18 


17 


79 


3.56 


1 


0.97 


0.95 






0.8464 


0.7453 


0.7310 






18.6 


13.4 


12.4 


65.4 


3.23 



25 

1 

1 

25 

30 

1 

1 

30 

24 

1 

1 

24 

17 

1 

1 

17 
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PT12 



23.5-128.6 



deep 


SSE 


4 


30 


30 


28 


20 


114 


4.27 


shoal 


k 

r 


1 


1 


1 


0.98 


0.963 






h = 325 


(k k 
s r 


1 


1 


0.8464 


0.7607 


0.7512 






a = 40 

r\ 


SSE 


4 


30 


25.4 


21.3 


15.0 


95.7 


3.91 



12 Aug 14 PT4 

24.6-126.7 



deep 


s 


1 


19 


31 


13 


8 


74 


3.44 


shoal 


k 

r 


1 


1 


0.975 


0.94 


0.91 






h = 266 


(k k 
s r 


1 
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58 


49 


39 
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39 


32 
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33 


25.3 
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18 


23 


17 
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81 
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44 


52 


36 
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30 
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18 
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1 
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1 

26 

26 
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1 

26 
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r 


1 


1 


1 
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5 
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1 


1 
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18 


12 


11 
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2.96 




a <10 
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36 
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00 Aug 01 
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